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Abstract 

We present a critical assessment of what can be learnt from the present data 
on inclusive polarized DIS. We examine critically some of the simplifying as- 
sumptions made in recent analyses and study in detail the question of the deter- 
mination of the gluon, strange sea and valence quark polarized densities. 

We have also carried out a new NLO QCD analysis of the world data. We 
find an excellent fit to the data and present our results for the polarized parton 
densities. 



1. Introduction. 



Deep inelastic scattering (DIS) of leptons on nucleons has remained the prime source 
of our understanding of the internal partonic structure of the nucleon and one of the 
key areas for the testing of perturbative QCD. Decades of experiments on unpolarized 
targets have led to a rather precise determination of the unpolarized parton densities. 
Spurred on by the famous EMC experiment |IJ at CERN in 1988, there has been a huge 
growth of interest in polarized DIS experiments which yield more refined information 
about the partonic structure. Many experiments have been carried out at SLAC @-@ 
and CERN 0-[lO] on proton, deuterium and 3 He targets, and there are major pro- 
grammes under way at SLAC (E155), DESY (HERMES 111) and CERN (COMPASS). 



In addition to the unpolarized structure functions Fi(x, Q 2 ) and F 2 (x,Q 2 ) there 
are two independent spin structure functions gi(x,Q 2 ) and g2(x,Q 2 ) and their 
unambiguous determination requires measurement of both the longitudinal asymmetry 
A\\ and the transverse asymmetry A± obtained with a target polarized parallel or 
perpendicular to the lepton beam direction, respectively. In the recent years there 
has been a great improvement in the quality of the data on the structure function 
g±(x, Q 2 ) , obtained from measurements using a longitudinally polarized target, and 
a big extension in the kinematic range x and Q 2 covered. Moreover it has become 
possible to present data in bins of (x, Q 2 ) rather than simply averaged over Q 2 at 
each x. The spin dependent structure function g2{x, Q 2 ) has now also been extracted 
|2], |5|, [12| from the data although with limited statistical precision compared to the g\ 
determination. 

The data at very small x has taught us that extrapolation of the measured values 
to x = is a subtle matter, so that the moments of the structure functions should 
not be considered as genuinely experimental quantities. 



Experiments on unpolarized DIS provide information on the unpolarized quark den- 
sities q(x, Q 2 ) and gluon density G(x, Q 2 ) inside a nucleon. Measurements of gi(x, Q 2 ) 
give us more detailed information, namely the number densities of quarks q(x, Q 2 )± 
and gluons G(x, Q 2 )± whose helicity is respectively along or opposite to the helicity 
of the parent nucleon. The usual densities are 

q(x, Q 2 ) = q + (x, Q 2 ) + q_(x, Q 2 ) , G(x, Q 2 ) = G+(x, Q 2 ) + <?_(*, Q 2 ) (1) 

and the new information is then contained in the polarized structure function g\(x, Q 2 ) 
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which is expressed in terms of the polarized parton densities 



Ag(x,Q 2 ) = q+ (x,Q 2 )-q4x,Q 2 ) , AG(x, Q 2 ) = G + (x, Q 2 ) - G-(x, Q 2 ) . (2) 
Several theoretical analyses [O - [ET1] based on the NLO perturbative QCD cal- 



culations [22] have sought to pin down the polarized parton densities. Each of them 



utilized the different data sets available at the time the analyses were performed. Only 



in the analyses [16. 19. 20 1 is essentially all the present data used, the exception being 
the very recent final E143 results ||). And none have completely used the informa- 
tion contained in the more detailed binning of the data in (x, Q 2 ). Moreover, there 
are differences in the assumptions to aid the analysis, differences in the choice of the 
renormalization scheme and differences in the form of the input parton densities and 
the value Qq at which they are determined. And finally, there are still significant dis- 
agreements about the results. To quote one example, Altarelli et al. [IfJ have obtained 



a significant polarized gluon density AG{x, Q 2 ) , whereas Bourrely et al. [21] claim 
that a perfectly acceptable fit can be obtained with AG(x, Q 2 ) = . 

Our aim in this paper is twofold: 

i) We discuss what can be learnt from the data, in theory and in practice, and 
examine the role played by the various assumptions used in the theoretical analyses. 

ii) We carry out a new study of the world polarized data. In addition to the data 
used in our previous analysis |TJJ the more accurate SLAC/E154 neutron data the 
new more precise SMC proton data |TUj which do not indicate a rise of gi at small 



x, the HERMES data [|TTJ and the final data results [Q of the E143 Collaboration at 
SLAC are now included. In trying to extract as well as possible the polarized parton 
densities we pay special attention to the observed scaling violations in the A\ data and 
use all the information contained in the more detailed binning of them in (x, Q 2 ). As in 
our previous work we use the following parametrization for the input polarized parton 
densities: 

Aq{x,Ql)=f(x)q{x,Q 2 ) , (3) 

in which we now utilize the new MRST set of unpolarized densities q(x, Qq) p3| . 
These parton densities account for the new, more precise HI and ZEUS deep inelastic 
scattering data, for the re-analysis of the CCFR neutrino data, for the inclusive prompt 
photon and large Et jet production in proton-proton collisions and for the charge asym- 
metry in Drell-Yan reactions. The MRST analysis leads to a value a s (M 2 ) = 0.1175 
in excellent agreement with the world average a s (M 2 ) = 0.118 ± 0.005 [f24j| . 
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2. Ambiguities and subtleties in determining the polarized parton den- 
sities 

There are several difficulties, specific to the polarized case, which make it much 
harder than in the unpolarized case to obtain reliable and unambiguous information 
about the polarized densities. 

2.1. What can be deduced in principle 

In the unpolarized case the separation of the contributions from partons of different 
flavours relies heavily upon the existence of both charged current neutrino and neutral 
current electromagnetic data. At present, and for some years to come, the information 
from polarized inelastic measurements will be limited to neutral current data. This 
raises the interesting question as to what one can hope to learn, both in theory and in 
practice. 

We have available data on g[ p \x,Q 2 ) and g^\x,Q 2 ) (or g± (x, Q 2 ) ) structure 
functions expressed as linear combinations of either the individual parton densities 

Au + Au, Ad + Ad, As + As (4) 

and AG or, equivalently, the 577(3) flavour combinations 

Ag 3 = (Am + Au) - (Ad + Ad) , 
Ag 8 = (Au + Au) + (Ad + Ad) - 2 (As + As) , 

AS = (An + Au) + (Ad + Ad) + (As + As) (5) 

and AG. (Note that Aq 8 is defined in such a way that its first moment is \/3 times the 
expectation value of the eighth component of the Cabibbo axial vector SU(3) current.) 

We are trying therefore to obtain information about four functions of x and Q 2 
on the basis of experimental data on the two independent functions g[ P \x,Q 2 ) 
and g[ n \x,Q 2 ). It is simpler to discuss the situation in terms of the contributions 
Ag 3 , Ag 8 , AS and AG. 

We have 

g p 1 {n) (x, Q 2 ) = \(^){5C NS [±^Ag 3 + ^Ag 8 ] + 5C S ® AS + 5C G ® AG} , (6) 

where 5Cns, 8Cs an d SCq are the non-singlet, singlet and gluon Wilson coefficient 
functions, respectively. In (||) ® denotes convolution with respect to x and 2/9 is the 
average value of the charge squared < e 2 > when the number of flavours Nf = 3. 
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Note that taking into account more than 3 active flavours does not change the main 
conclusions in this and the next section. 



Let us firstly suppose that we have perfect data for a range of x and Q 2 , and that 
we try to determine the functions Aq 3 , Aq s , AS and AG. Then Aq 3 is determined 
uniquely and trivially since 

g{{x, Q 2 ) - g?(x, Q 2 ) = ^5C NS ® Ag 3 (7) 
We are then left with 

g{{x, Q 2 ) + g?(x, Q 2 ) = \[\sC NS ® Ag 8 + 6C S ® AS + 6C G ® AG] . (8) 

It is the difference in the Q 2 -evolution of the three terms on the RHS of (H) that 
enables them to be determined separately. Indeed, by studying the first and higher 
derivatives of the LHS of @ with respect to Q 2 at Q 2 = Ql, and using the evolution 
equations, one can prove that Aq 8 , AS and AG are uniquely determined at Q 2 = Ql- 

It is immediately clear, given the limited range of Q 2 available and the fact that 
the data are not perfect and have errors, that the separation of Aqs, AS and AG 
from each other will not be very clearcut. Nonetheless, in principle, the data fix 
Ag 3 , Ag 8 , AS and AG or, equivalently, via (f|), Au + Au, Ad + Ad, As + As = 2As 
and AG. But any hope of a successful analysis will depend upon finding simple enough 
parametrizations of these quantities at Q 2 = Q\- 



2.2. Valence and sea 



It is clear from the above that the inclusive (electromagnetic current) data give 
no information about the valence parts Aq v of the quark densities. It is of interest, 
however, to know the Ag^ and Aq or equivalently, Ag and Aq , since they 
play a distinctive role in other types of experiment, e.g. polarized semi-inclusive DIS, 
polarized Drell-Yan reactions etc. Indeed, an attempt has been made to extract the 
polarized valence densities from the semi- inclusive data [P5| , but the quality of the 
present data precludes an accurate determination of these densities. It is therefore 



important to make a combined analysis |zb] of both the semi-inclusive and inclusive 
DIS data. 

Further, for the unpolarized densities simple parametrizations are normally given 
for the valence and sea quark densities. One reason for this is that one expects q v and 
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q to have simple, but very different, behaviours as x — > 0, and this feature is lost when 
dealing with q or q + q. 

There is some point, therefore, in wanting to deal directly with the valence and sea 
densities and to this end it has been common practice to make some model assumptions 
about the polarized sea jl3[ [14J] and []17| |- [ [2l"|l , which then allows a determination of 



the valence parts. For example, the apparently innocuous assumption of a flavour 
independent polarized sea 

As = Au = Ad (9) 

implies that Au v and Ad v are determined, since the data fix As, Au v + 2Au and 
Ad v + 2 Ad. 

We believe it is important to study the consequences of assumptions like Eq. 
Consider therefore the family of assumptions at Q 2 = Ql 

Au = Ad = AAs , (10) 

where A is a parameter. 

Given that the data fix Aq^g, AS and AG and that 

As = 1(AE-Ag 8 ) , (11) 
o 

we see that the result for As should not change as A is varied. This provides a serious 
test for the stability of the analysis. 

Further the dependence of the valence densities upon A is given by 

Au v = ^[Aq 3 + Ag 8 - 4(A - l)As] , 

Ad v = ^[-Ag 3 + Ag 8 -4(A-l)As] , (12) 

so that they are sensitive to the assumption about the sea. On the other hand, if the 
analysis is correct, neither Aq 3 8 nor AS(As) nor AG should change as A is varied. 
It should be noticed that Eqs. ( PH) and flT2"|) are only valid at Q 2 = Ql- The equality 



(10) and therefore Eq. (p!2|) will be (marginally) broken because of the different NLO 
evolution of the different sea quarks for Q 2 > Q^. 



It has sometimes been claimed in the literature |19| that as a result of the fits the 
sea always turns out to be flavour symmetric. Clearly, from the above, one can learn 
absolutely nothing from the data about A-u and Ad and any \ 2 dependence on A 
must be an artifact of the fitting procedure. 
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2.3. Simplifying assumptions 



In order to limit the number of parameters being fitted, it is always necessary 
to take simple functional forms for the parton densities, and sometimes additional 
simplifying assumptions are made. We have already discussed in the previous section 



the assumptions regarding the sea quarks (see Eqs. @ and (|10D ) widely used in the 
literature. Another somewhat arbitrary one (for its more recent use see Altarelli et al. 
HID!), is to assume that 

Aq 3 (x,Q 2 ) = CAq s (x,Q 2 ) , (13) 

where C is a constant. 

This is, of course, perfectly compatible with the evolution equations, decreases sub- 
stantionally the number of parameters being fitted, but has no physical justification at 
all. This assumption leads to a better determination of the rest of the independent par- 
ton densities, but it is not clear whether the values and errors of the parameters are not 



thereby distorted. It is important to check to what extent ( |i3|) is compatible with the 
results of other theoretical analyses of the data in which this approximation is not used. 



2.4- Scheme dependence 



It is well known that at NLO and beyond, the parton densities become dependent 
upon the renormalization (or factorization) scheme. In the unpolarized case the most 
commonly used are the MS, MS and DIS schemes and parton densities in different 
schemes differ from each other by terms of order a s (Q 2 ) , which goes to zero as Q 2 
increases. 

There are two significant differences in the polarized case: 

i) The singlet densities AE(x, Q 2 ), in two different schemes, will differ by terms of 
order 

a s (Q 2 )AG(x,Q 2 ) , (14) 

which appear to be of order a s . But it is known ]37L ES| that, as a consequence of the 
axial anomaly, the first moment 

f 1 dxAG(x,Q 2 ) oc [a s (Q 2 )Y 1 , (15) 
Jo 

grows in such a way with Q 2 as to compensate for the factor a s (Q 2 ) in fll4|) . Thus 
the difference between AE in different schemes is only apparently of order a s (Q 2 ), and 
could be quite large. 
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ii) Because of ambiguities in handling the renormalization of operators involving 
75 in n dimensions, the specification MS does not define a unique scheme. Really 



there is a family of MS schemes which, strictly, should carry sub-label indicating 



how 75 is handled. What is now conventionally called MS is in fact the scheme due 
to Vogelsang and Mertig and van Neervenen f22|, in which the first moment of the 



non-singlet densities is conserved, i.e. is independent of Q 2 , corresponding to the 
conservation of the non-singlet axial-vector Cabibbo currents. 

Although mathematically correct it is a peculiarity of this factorization scheme 
that certain soft contributions are included in the Wilson coefficient functions, rather 
than being absorbed completely into the parton densities. As a consequence, the first 
moment of AS is not conserved so that it is difficult to know how to compare the DIS 
results on AE with the results from constituent quark models at low Q 2 . 



To avoid these idiosyncrasies Ball, Forte and Ridolfi introduced what they 
called the AB scheme, which involves a minimal modifications of the MS scheme, and 
for which 

AZ(x,Q 2 ) AB = AE(x, Q 2 )ms + Nf^Q-t f 1 ^-AG(y, Q 2 )^ , 

ZTT Jx y 

AG(x,Q 2 ) AB = AG(x,Q 2 )j^ (16) 
or, in the Mellin n-moment space, 



AZ(n,Q 2 ) AB = AZ(n,Q 2 )j^+ Nf^^AG^Q 2 )^ , 
AG(n, Q 2 )ab = AG(n, Q 2 )jjs ■ (17) 



That AS(n = 1) AB is independent of Q 2 to all orders follows from the Adler- 
Bardeen theorem [^J. 

The singlet part of the first moment of the structure function g 1 

rj s) (Q 2 )= f'dxg^^Q 2 ) (18) 
Jo 

then depends on AE and AG only in the combination 

a (Q 2 ) = AS(1, Q 2 )j^ = AE(1) AB - iV^AGfl, Q 2 ) (19) 
and the unexpectedly small value for the axial charge found by the EMC [l|] , which 



triggered the "spin crisis in the parton model" |30|], can be nicely explained as due 
to a cancellation between a reasonably sized AS(1) and the gluon contribution. Of 
importance for such an explanation are both the positive sign and the large value (of 
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order 0(1)) for the first moment of the polarized gluon density AG(1,Q 2 ) at small 
Q 2 ~ 1 - 10 GeV 2 . Note that what follows from QCD is that |AG(1,Q 2 )| grows with 
Q 2 (see Eq. (^)) but its value at small Q 2 is unknown in the theory at present and 
has to be determined from experiment. 



Although the AB scheme corrects the most glaring weakness of the MS scheme, it 
does not consistently put all hard effects into the coefficient functions. As pointed out 
in one can define a family of schemes labelled by a parameter a: 

where 

z qG (x; a) = N f [(2x - l)(a - 1) + 2(1 - x)] , (21) 

in all of which (|19|) holds, but which differ in their expression for the higher moments. 
(The AB scheme corresponds to taking a = 2). 

Amongst these we believe there are compelling reasons to choose what we shall call 
the JET scheme (a = 1), i.e. 

z J J T = 2N f (l - x) . (22) 

This is the scheme originally suggested by Carlitz, Collins and Mueller and also 
advocated by Anselmino, Efremov and Leader ||32|| .' In it all hard effects are absorbed 
into the coefficient functions. In this scheme the gluon coefficient function is exactly 
the one that would appear in the cross section for 

W -> jet(k T ) + jet{-k T ) + X , (23) 

i.e., the production of two jets with large transverse momentum ky and — k^, respec- 
tively. 

More recently Miiller and Teryaev |33| have advanced rigorous and compelling ar- 
guments, based upon a generalization of the axial anomaly to bilocal operators, that 
removal of all anomaly effects from the quark densities leads to the JET scheme. Also 
a different argument by Cheng |34| leads to the same conclusion. (Cheng calls the JET 
scheme a chirally invariant (CI) scheme.) 



The transformation from the MS scheme of Mertig, van Neerven and Vogeslang to 
the JET scheme is given in moment space by 

A£(n,Q 2 ) J£T = AZ(n J Q 2 ) J ^ + 2N f 2 2^ 1) AG(n,Q 2 ) JT s , 

tThere is misprint in Eq. (8.2.6) of ||. The term ln( l ~ x J? ' ) should be [In ^^Jj' ) - 1]. 
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AG(n,Q 2 ) JET = AG^Q 2 )^. (24) 
Of course, ([T7|) and (PI ) become the same for n = 1. 



In this paper we carry out the fitting procedure in the MS scheme and the results 



can then be transformed to the other schemes via (|1^) and ([24]) • However it will be 
important to carry out the fitting in the other schemes as a check on the stability of 
the whole analysis |35 . 



3. Method of analysis and input parton distributions 

The spin dependent structure function of interest, g^(x, Q 2 ) , is a linear combina- 
tion of the asymmetries AW and (or the related virtual photon-nucleon asymmetries 
Ai 2 ) measured with the target polarized longitudinally or perpendicular to the lepton 
beam, respectively. Neglecting as usual the subdominant contributions (see for exam- 
ple M), Ai(x, Q 2 ) can be expressed via the polarized structure function g±(x,Q 2 ) 



as 



where 



A-(x, Q 2 ) ^ (1 + 1 2 )§W^ = fplfyt 2 ^ 1 + RN ^ <? 2 )] > ( 25 ) 
R N + 1 = (1 + 7 2 )F 2 7V /2xFf (26) 



and F± and FJ? are the unpolarized structure functions. In ( 25|) the kinematic 
factor 7 2 is given by 



7" = ^. (27) 



AM 2 N x 2 

It should be noted that in the SLAC kinematic region 7 cannot be neglected. 

In some cases the theoretical analyses of the data are presented in terms of gi(x, Q 2 ) 
as extracted from the measured values of A±(x, Q 2 ) according to (p5|), using various 
parametrizations of the experimental data for F2 and R. 

As in our previous analysis we follow the approach first used in jTJ]], in which the 
next-to-leading (NLO) QCD predictions for the spin-asymmetry A^(x, Q 2 ) are con- 
fronted with the data on A^a;,^ 2 ) , rather than with the gi(x,Q 2 ) derived by 
the procedure mentioned above. The choice of A± should minimize the higher twist 
contributions which are expected to partly cancel in the ratio (|25|), allowing use of data 
at lower Q 2 . Bearing in mind that in polarized DIS most of the small x data points 
are at low Q 2 , a lower than usual cut is needed ( Q 2 > 1 GeV 2 ) in order to have 
enough data for the theoretical analysis. We believe that in this approach such a low 
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Q 2 -cut is more justified. 



In NLO approximation 



g 1 (x,Q ) NLO ,, 01 

A x {x, Q ) NLO = (1+7) pN , n2 , • 

r 1 [x,LJ )nlo 

In (£D N = p, n and d = (p + n)/2 . 

To calculate A±(x, Q 2 )nlo in NLO QCD and then fit the data we follow the 
same procedure described in detail in our paper Here we will recall only the main 
points. 

The Q 2 evolution, in NLO QCD approximation, is carried out for the n-space 
moments of the polarized quark and gluon densities. Then using the known NLO 
expressions for the moments of the Wilson coefficients 5C q (n,a s ) and SCcin, a s ) 
(see e.g. |L4|]) one can calculate the moments of the structure function 



i N f 

M N (n,Q*) = ~Y,e 2 q 
/ q 



5C q (n)(Aq(n, Q 2 ) + Aq(n, Q 2 )) + -^-8C G (n)AG(n, Q' 

iv f 



(29) 

As already mentioned above, all calculations are performed in the MS scheme. To 



account for heavy quark contributions we use the so-called fixed-flavour scheme |36], [14 
and set the number of active flavours in (p9|) Nf = 3. In contrast to our previous 



analysis [0, we now use for the values of the QCD parameter Aj^g : Aj^g(nf = 3) = 
353 MeV and A M (n/ = 4) = 300 MeV , which correspond to a s (M 2 ) = 0.1175 , 
as obtained by the MRST analysis |23] of the world unpolarized data, in excellent 
agreement with the world average a s (M 2 ) = 0.118 ± 0.005 [[24] . 



Finally, to reconstruct the spin structure functions gi(x,Q 2 ) in Bjorken x-space 



from their moments (|29|) with the required accuracy, we use the Jacobi reconstruction 
method [[37], |38|. Note that in this method the structure functions are given analytically. 

The same procedure has been used to calculate the unpolarized structure functions 
Ff(x, Q 2 )nlo from their moments. 

We choose the input polarized densities at Qq — 1 GeV 2 in the form: 

xAu v {x,Qq) = Tj u A u x u xu v {x^Q^) , 

xAd v (x,Ql) = r] d A d x ad xd v (x,Ql) , 

xASea(x, Ql) = i]sAsx as xSea(x, Ql) , 

xAG(x,Ql) = Vg A g x a *xG(x,Q 2 ) (30) 
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where on R.H.S. of (|30|) we have used the MRST unpolarized densities |23 . 

Guiding arguments for such an ansatz are simplicity (not too many free parameters) 
and the expectation that polarized and unpolarized densities have similar behaviour 



at large x. In ( pOf ) the parameters a/ account for the difference of the low- a; behaviour 
between the polarized and unpolarized parton densities. The normalization factors 
Af are determined in such a way as to ensure that the first moments of the polarized 
densities are given by r]f. 

In the previous section we explained why we chose to deal with valence and sea 
quarks instead of their singlet and non-singlet combinations. For the polarized light 
and strange sea quark densities at Qq = 1 GeV 2 we adopt the assumption ([Top. Then 

a- -a- ASea 

AS = Aq = 2(2XTT) ' * = 2(2ATl)' (31) 

where r] s is the first moment of the strange sea parton density As. 

We would like to emphasize once more that in contrast to the valence quark densi- 
ties, As should not depend on the flavour sea decomposition, i.e. on A in our case (see 
Eq. (|TTD), and as will be seen below, our numerical results confirm this. 

The first moments of the valence quark densities r] u and % are fixed by the octet 
hyperon (3 decay constants |5U 



g A = F + D = 1.2573 ± 0.0028, a s = 3F — D = 0.579 ± 0.025 . (32) 

to be 

Vu = 0.918 -2(A-l)As(l,Qg), % = -0.339 - 2(A - l)As(l, Qg). (33) 

In the case of SU(3) flavour symmetry of the sea (A = 1) we take 

r] u = 0.918 , r] d = -0.339 . (34) 

The rest of the parameters in (pQj) 

{a u , a d , T]s, a s ,% ,a g } , (35) 

have to be determined from the best fit to the Af(x, Q 2 ) data. 

In some papers |LJ| ^1|] gA, and in others (see e.g. [0]) both g A and ctg have been 
taken to be free parameters determined by the best fit to the inclusive polarized DIS 
data. We do not favour such an approach because the values of these quantities, espe- 
cially gA, are determined from much more precise experiments and using them improves 
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the accuracy with which we can determine the polarized densities. 



Results of Analysis 

In this section we present the results of our fits to the present experimental data 
on A?(x,Q 2 ) : EMC proton data SLAC E142 neutron data 0, SLAC E154 



neutron data ||, SMC combined proton data [10|, the SMC deuteron data || which 
are combined data from the 1992, 1994 [[| and 1995 runs, HERMES neutron data JTTJ] 
and the final SLAC E143 results [| on g{/Ff and gf/Ff . The data used (354 
experimental points) cover the following kinematic region: 

0.004 < x < 0.75, 1 < Q 2 < 72 GeV 2 . (36) 

As mentioned in the Introduction, in contrast to the other analyses, we fit all pos- 
sible (x, Q 2 ) data rather than ones averaged over Q 2 within each x-bin. We denote this 
combined fit to the (x, Q 2 ) data on presented by E142, E143 and SMC collabora- 
tions and the averaged A± data given by EMC, E154 and HERMES as Fit A. Since for 
most of these data (E142, SMC) the systematic errors are not published, in Fit A only 
statistical errors are taken into account. The results of the fit to the averaged A^ data 
alone (118 experimental data points) given by all the collaborations mentioned above 
are also presented (Fit B). In this case the total (statistical and systematic) errors are 
included in the analysis. "Higher twist" corrections are not included in the present 
study. As already discussed above, in the approach used their effect is expected to be 
negligible. 



The numerical results of the fits (A = 1) are listed in Table 1. Note that r/s = rjs/Q in 
the SU(3) symmetry case. The dependence of the results on the flavour decomposition 
of the sea will be discussed below in detail. 

It follows from our analysis that the value of a g can not be well determined, i.e. 
the existing data do not constrain the behaviour of the polarized gluon density at 
small x. For that reason the fits to the data were performed at different fixed val- 
ues of a g in the range:0 < a g < 1. In Fit A the change of x 2 /DOF value from 
X 2 /DOF{a g = 0) = 0.917 to x 2 /DOF{a g = 1) = 0.910 is negligible. The same 
conclusion is true for the fits to the average A\ data. In Table 1 we present the results 
of the fits corresponding to a g = 0.6. 

In Fig. 1 the SLAC/E143 and SMC data on A\ and A\ vs Q 2 for different x-bins 
are compared to our best Fit A. The NLO results for the averaged asymmetries A± 
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(Fit B) are shown in Fig. 2. It is seen from the values of x 2 /DOF and Figs. 1 and 
2 that the NLO QCD predictions are in a very good agreement with the presently 
available data on , as well as with the corresponding g^(x, Q 2 ) data (see Figs. 6a 
and b). We would like to draw special attention to the excellent fit to the E154 neutron 
data (see Fig. 2c), the most accurate polarized DIS data at present (x 2 = 1-6 for 11 
experimental data points). 



Table 1. Results of the NLO QCD fits to the world A? data (Q 2 = 1 GeV 2 ). 
For Fit A errors are statistical, for Fit B, total. a g = 0.6 (fixed). 



Parameters 


Fit A 


Fit B 


DOF 


354 - 5 


118 - 5 


x 2 


318.4 


86.1 


X 2 /DOF 


0.912 


0.762 


a u 


0.250 ± 0.023 


0.255 ± 0.028 


ad 


0.231 ± 0.088 


0.148 ± 0.113 


a s 


0.576 ± 0.152 


0.817 ± 0.223 


Vs 


- 0.054 ± 0.012 


- 0.049 ± 0.005 


% 


0.34 ± 0.24 


0.82 ± 0.32 




0.253 ± 0.079 


0.287 ± 0.041 


AS(1)ab 


0.332 ± 0.096 


0.476 ± 0.084 



One can see from Fig. 1 that the accuracy and the presently measured kinematic 
region of the data do not allow a definite conclusion about the scaling violations in 
Ai(x, Q 2 )- It is obvious that more precise data and an extension of the measured 
range to smaller x and larger Q 2 are needed to answer the question about the Q 2 de- 
pendence of virtual photon spin asymmetry Af. 

The extracted valence, strange and gluon polarized distributions at Q 2 = 1 GeV 2 
are shown in Fig. 3a: Fit A (solid curves) and Fit B (dashed curves). Except for 
the gluon these two sets of densities coincide almost exactly in the measured x-region. 
The polarized gluon densities extracted in Fit A and Fit B are a good illustration of 
how large the uncertainty is in determining the gluon density from the present data. 
Although the values of their first moments i] g are in agreement within two standard 
deviations, their central values differ by a factor of more than two. (Note that in Fit 
A and B the same form of the initial parton densities has been used.) 
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In Fig. 3b and Fig. 4 the nonsinglet combinations x Aq% , xAq$ and their ratio 
Aq 3 /Aq 8 are shown, respectively. The ratio Aq 3 /Aq$ is shown at Q 2 = 1 GeV 2 , but it 
does not change with Q 2 because the same evolution holds for both xAq 3 and xAq$. 
We find a significant deviation from the approximation fllTf ) used in jL6| by Altarelli 
et al. Note that the Q 2 evolution of the nonsinglets is the same in all the schemes 
discussed in Section 2.4. 



Let us now examine how the assumption ( |i~0| ) about the flavour decomposition of 
the sea influences our results. As pointed out in Section 2.2 the strange quark density, 
and in particular, its first moment As(1,Qq) should not change as A is varied. The 
results of the fits to the averaged data (Fit B) using different values of A are 
presented in Table 2. 

Table 2. The results for the first moments of the polarized distributions at 
Q 2 = 1 GeV 2 using the assumption ([K]) about the flavour decomposition of 
the sea. 



A 


x 2 


- Vs 


% 


Vu 


- Vd 


0.5 


85.97 


0.049 ± 0.004 


0.78 ± 0.31 


0.869 ± 0.013 


0.388 ± 0.013 


1.0 


86.11 


0.049 ± 0.005 


0.82 ± 0.32 


0.918 ± 0.013 


0.339 ± 0.013 


2.0 


86.08 


0.050 ± 0.007 


0.86 ± 0.34 


1.018 ± 0.019 


0.239 ± 0.019 


3.0 


86.12 


0.048 ± 0.004 


0.93 ± 0.35 


1.110 ± 0.020 


0.147 ± 0.020 



It is clear from the table that \ 2 an d the central values of r] s = As(l, Qo) and 7] g = 
AG(1, Ql) practically do not change as A varies. We regard this fact as a very good test 
of the stability of our analysis. We thus conclude, somewhat in disagreement with Ref. 



16| that the separation into valence and sea contributions need not introduce biases 



into the fit provided sufficient care is taken. Of course, we have also demonstrated very 
clearly that the Au v and Ad v valence quark densities are sensitive to the assumptions 
about the sea. This dependence is shown in Fig. 5. 

What follows from our analysis is that one can use for input distributions a parametriza- 
tion in terms of valence and sea quarks. With a correct fitting procedure the strange 



sea (or singlet, see (|Tl"D) and gluon distributions defined from the inclusive data do 
not depend on the assumption about the sea and therefore, one can use them to test 
the remarkable relation fll9"|). On the contrary, as emphasized in Section 2.2, elec- 
tromagnetic DIS does not fix the valence quark densities, which are sensitive to the 
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assumption about the sea. This implies that is somewhat meaningless to claim |26| 
that "the present semi-inclusive data alone fail to define a Ad v distribution consistent 
with those extracted from inclusive data" . 



In Table 1 we also present our results for the first moments of the polarized gluon 



and singlet quark densities in the MS (determined from the fit) and AB (calculated 
by Eq. flT7D) schemes. The value of r\ g (Fit B) 



rj g = AG(1, 1 GeV 2 



'MS 



AG(1, 1 GeV 2 ) A B = 0.82 ± 0.32 



(37) 



is in a good agreement with the one obtained in Fit A of Ref . |16|] : 

AG(1,1 GeV 2 ) A B = 0.95 ±0.18 (38) 

using almost the same data. 
Our result for AE(1)ab 

AE(1)ab = 0.476 ± 0.084 (39) 
is in agreement within two standard deviations with its constituent value 0.6 |30] and, 



within errors, coincides with the value determined in |]I6|| 



AE(1) AB = 0.405 ±0.032 . 



(40) 



Finally, let us turn to the first moments T^(Q 2 ) of the spin structure function 
g^(x,Q 2 ). Using our results of the fits for the input polarized parton densities these 
quantities have been calculated for different values of Q 2 using Eq. (|29| ) for n = 1. The 
results are presented in Table 3. 

Table 3. Determination of the first moments T^(Q 2 ) of gi(x,Q 2 ) 
using the results of our NLO QCD fits to the world data. 





Fit A 


Fit B 


r?(5 GeV 2 ) 
17(5 GeV 2 ) 
rf(5 GeV 2 ) 
r?(10 GeV 2 ) 
17(10 GeV 2 ) 
17(10 GeV 2 ) 


0.134 

- 0.056 
0.036 
0.136 

- 0.057 
0.036 


0.138 

- 0.053 
0.039 
0.139 

- 0.054 
0.039 
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The values of the first moments of the polarized parton densities and therefore, the 
corresponding values of the first moments of gf are very sensitive to the assumed small 
x-behavior of the input parton densities (see e.g. || fLOp. Parton densities which give 
the same results for the physical quantities in the measured x region can lead to rather 
different behaviour at very small x. The dependence of the the moments of the physical 
quantities in polarized DIS on different assumptions about the small x behaviour of 
the input parton densities has been studied in detail in fT6fl . 

The simple Regge behaviour of the unpolarized and polarized structure functions as 
x — > 0, which was for many years a guiding principle in our choice of the starting ansatz 
for the parton densities, differs from that predicted in QCD [[y], |42|. Moreover, the 
unpolarized DIS experiments at HERA have shown that the simple Regge extrapolation 
of the structure functions at small x is not valid at large Q 2 . There is an indication || 
that this is also true in the polarized case. The question of the small x behaviour in 
the polarized case remains an open one and is a serious challenge to both experiment 
and theory. 

For illustration we compare our predictions for 17(10 GeV 2 ) and 17(5 GeV 2 ) 

r?(10 GeV 2 ) = {0.136 (Fit A), 0.139 (Fit B)} , 
17(5 GeV 2 ) = {-0.057 (Fit A), -0.054 (Fit B)} (41) 



with their "experimental" values |T0|] and |T7] 



17(10 GeV 2 ) = 0.130 ± 0.006(stat) ± O.OOS(syst) ± O.OU(evol) , (42) 

17(5 GeV 2 ) = -0.058 ± 0.004(stat) ± 0.007(syst) ± 0.007(evol) , (43) 

the latter chosen because the moments have been determined by extrapolating the 
polarized structure functions into the unmeasured x region using QCD. It is seen from 
Eqs. (|4T| ) - (|43|) that our predictions for 17 and T™ are in a good agreement with their 
"experimental" values. 



5. Conclusion 

We have performed a next-to leading order QCD analysis (MS scheme) of the 
world data on inclusive polarized deep inelastic lepton-nucleon scattering. The QCD 
predictions have been confronted with the data on the virtual photon-nucleon asym- 
metry A^(x, Q 2 ) , rather than with the polarized structure function g±(x,Q 2 ) , in 
order to minimize the higher twist effects. In this paper, for the first time, is utilized 
the full world data set on A± with its detailed binning in (x,Q 2 ). Using the simple 
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parametrization (]30|) (with only 5 free parameters) for the input polarized parton den- 
sities it was demonstrated that the polarized DIS data are in an excellent agreement 
with the pQCD predictions for Ai(x, Q 2 ) and the spin-dependent structure function 
gi(x, Q 2 ) . However, the accuracy and the presently measured kinematic region of the 
data do not allow a definite conclusion about scaling violations in A^(x, Q 2 ). 

We have studied the consequences of different simplifying assumptions, usually 
made in recent analyses to aid the extraction of the polarized parton densities from the 
data. It was shown that whereas the valence quark densities determined from inclusive 
polarized DIS data, are sensitive to the assumptions about the flavour decomposition 
of the sea, the extracted strange sea and gluon densities in our analysis do not depend 
on such an assumption and can therefore be used to test the remarkable relation ([TDD. 



We have found a significant deviation from the approximation ( ]T3| ) used in Ref. [IB 
that the nonsinglet quark distributions Aq 3 (x,Q 2 ) and Aq 8 (x,Q 2 ) have the same 
shape at fixed Q 2 . 

Although the quality of the data has significantly improved via the recent experi- 
ments of the SLAC/E154 Collaboration and the final SMC results on Af, the uncer- 
tainty in determining the polarized gluon density is still very large. 

Despite the great progress of the past few years it is clear that in order to test 
precisely the spin properties of QCD, more accurate inclusive DIS polarized data and 
an extension of the measured range to smaller x and larger Q 2 are needed. We 
hope the current DIS experiments at HERA will help in clarifying the situation. In 
addition, semi-inclusive and charged current data will be very important for a precise 
determination of the polarized parton densities and especially, for an accurate flavour 
decomposition of the polarized quark sea. There is some progress in this direction . 



Finally, a direct measurement of AG(x, Q 2 ) in processes such as J/ip production in 
lepton-hadron scattering with a polarized beam will answer the important question 
about the magnitude of the first moment of the gluon density AG(x, Q 2 ). 
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Figure Captions 

Fig. la-lb. Ai and Af vs Q 2 for different x-bins for the proton E143 [f§ and SMC 
data [1C] and for deuteron E143 |3J and SMC data ||. Only statistical errors are shown. 
The solid curves correspond to Fit A described in the text. The E143 data on gi/Fi 
are multiplied by the kinematic factor (1 + j 2 ) . 



Fig. 2a-2c. Comparison of our NLO results (Fit B) for averag ed A?(x,Q 2 ) with the 
experimental data at the measured x and Q 2 values. Errors bars represent the total 
error. 



Fig. 3a-3b. Next-to-leading order input polarized parton distributions at Q 2 = 
1 GeV 2 (A = 1). Solid and dashed curves in Fig. 3a correspond to Fit A and Fit B, 
respectively. 

Fig. 4. Comparison between our result (solid curve) for the ratio Aq^(x) j Aqg(x) and 
approximation flT3] ) used in Ref. |TB[ (dashed line). 

Fig. 5. Polarized valence quark densities at Q 2 = 1 GeV 2 for different values of A 
(see Eq. ([T0|)). The solid curves correspond to an SU (3) flavour symmetric sea (A = 1). 

Fig. 6a-6b. Comparison of our NLO results for g^(x,Q 2 ) (6a) and xg^(x, Q 2 ) 
(6b) with SMC [g [10[ and SLAC/E154 data at the measured values of Q 2 . Error 
bars represent the total error. 
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